
** note only works if run after table 4. 

frame change default 





* frankenmodels JD plus one of predicted_jd predicted_dd_21 predicted_ks_21 
constraint 1 predicted_jd + predicted_honesty=1
constraint 2 predicted_jd + predicted_dd=1
constraint 3 predicted_jd + predicted_ks=1

gen predicted_dd = predicted_dd_21
gen predicted_ks = predicted_ks_21
est clear

* need to fix BIC and AIC to take account of number of parameters
eststo: cnsreg reports_21 predicted_jd predicted_honesty, constraint(1) nocons
local parameter_k 1
estadd scalar RMSE_raw  =  e(rmse) * sqrt(e(df_r) / e(N))
estadd scalar AIC = (2*`parameter_k')-2*e(ll)
estadd scalar BIC = `parameter_k'*ln(36)   - 2*e(ll)
eststo: cnsreg reports_21 predicted_jd predicted_dd  , constraint(2) nocons
estadd scalar RMSE_raw  =  e(rmse) * sqrt(e(df_r) / e(N))
local parameter_k 2
estadd scalar AIC = (2*`parameter_k')-2*e(ll)
estadd scalar BIC = `parameter_k'*ln(36)   - 2*e(ll)
eststo: cnsreg reports_21 predicted_jd predicted_ks , constraint(3) nocons
estadd scalar RMSE_raw  =  e(rmse) * sqrt(e(df_r) / e(N))
local parameter_k 3
estadd scalar AIC = (2*`parameter_k')-2*e(ll)
estadd scalar BIC = `parameter_k'*ln(36)   - 2*e(ll)

replace predicted_dd = predicted_dd_22
replace predicted_ks = predicted_ks_22
eststo: cnsreg reports_22 predicted_jd predicted_honesty, constraint(1) nocons
estadd scalar RMSE_raw  =  e(rmse) * sqrt(e(df_r) / e(N))
local parameter_k 1
estadd scalar AIC = (2*`parameter_k')-2*e(ll)
estadd scalar BIC = `parameter_k'*ln(36)   - 2*e(ll)
eststo: cnsreg reports_22 predicted_jd predicted_dd  , constraint(2) nocons
estadd scalar RMSE_raw  =  e(rmse) * sqrt(e(df_r) / e(N))
local parameter_k 2
estadd scalar AIC = (2*`parameter_k')-2*e(ll)
estadd scalar BIC = `parameter_k'*ln(36)   - 2*e(ll)
eststo: cnsreg reports_22 predicted_jd predicted_ks , constraint(3) nocons
estadd scalar RMSE_raw  =  e(rmse) * sqrt(e(df_r) / e(N))
local parameter_k 3
estadd scalar AIC = (2*`parameter_k')-2*e(ll)
estadd scalar BIC = `parameter_k'*ln(36)   - 2*e(ll)


esttab * using "output/table5.tex", stats(AIC BIC  RMSE_raw)  star(* .1 ** .05 *** .01) tex  replace